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We present a simple model in which dark matter couples to the standard model through a light 
scalar intermediary that is itself unstable. We find that this model has several notable features, and 
allows a natural explanation for a surplus of positrons and electrons, but no surplus of anti-protons, 
as has been suggested by early data from PAMELA and ATIC. Moreover, this model yields a very 
small nucleon coupling, well below the direct detection limits. In this paper we explore the effect of 
this model in both the early universe and in the galaxy. 



I. INTRODUCTION 

There is very strong evidence that nonbaryonic dark 
matter comprises about 21% of the present energy den- 
sity of the universe. Identification of the dark compo- 
nents poses one of the strongest challenges in particle 
physics and astrophysics, and will almost certainly re- 
quire an extension to the standard model. Indeed, there 
is a rich literature on dark matter candidates that arise 
more or less naturally from various extensions. 

A number of recent experiments have reported obser- 
vations that, if confirmed in their final analysis, may be 
the first strong non-gravitational hints to the identity of 
dark matter. The PAMELA collaboration has reported 
that their satellite-born experiment has observed an ex- 
cess in the ratio of positrons to the total expected electron 
plus positron background. The spectrum they reported 
departed from background at 10 GeV and rose monoton- 
ically to 100 GeV[T]. Recently the PAMELA experiment 
has provided an updated spectrum that incorporates ad- 
ditional data and a new statistical method [Ij. We will 
use the "beta-fit" from that paper. Their final data re- 
lease is expected to have the positron spectrum measured 
to about 300 GeV [HIS]. The balloon-based ATIC exper- 
iment measured the electron spectrum in the range 10 
GeV to 2 TeV. Their spectrum, given in [3], seems to 
show a significant excess from 300 GeV to 600 GeV, with 
a peak at roughly 500 GeV. This data is consistent with 
published data from another balloon experiment, PPB- 
BETS, which measured electrons up to about 1 TeV [5], 
and with the earlier hints of an excess at the HEAT ex- 
periment [^. The Fermi Space Telescope has released 
a combined electron and positron spectrum with excel- 
lent statistics across the ATIC band, from 20 GeV to 
1 TeV [7 . Their measurement does not agree with ex- 
cess found by ATIC. The spectrum does appear to have 
an excess of events with a peak near 400 GeV, though the 
authors emphasize that the spectrum may be fit with a 
simple power law. 

A number of groups have studied the consistency of 
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these experiments with various theories of dark matter. 
Some groups have considered the signals generated by 
specific models, including Minimal Dark Matter [Sj the 
two-component A'dmMSSM model [9 , supersymmetry 
[101 lll| and a model with a nonthermal production mech- 
anism [12]. Other groups have instead fit the data to 
general models as well as specific models [131 [T3] . Re- 
cently, Arkani-Hamed et. al. have proposed an interest- 
ing dark matter model that takes advantage of the Somer- 
field mechanism to achieve enhanced annihilation '151, 
following earlier work in refs. [14) . Their model, as well 
as the model of Pospelov et al. [THl [T7] , share some char- 
acteristics with the one discussed in this paper due to 
the inclusion of a light boson to generate the lepton spec- 
trum. Earlier work by Cholis et al. [TS] studied a scenario 
in which dark matter annihilates through a light inter- 
mediary. A study by March-Russel et al. [19] considers a 
supersymmetric dark matter model with a similar mech- 
anism of communication with the standard model. Some 
early work on which galactic signals are generated by par- 
ticular models appears in [20H22] . Many of these models 
associate the dark matter with a new sector containing 
multiple exotic new particles. 

In this paper we present a dark matter model that ex- 
plains the experimental observations in a simple way and 
differs significantly from the scenarios mentioned above. 
We assume that there is a new sector of particles that 
interact weakly with standard model particles, but may 
have stronger interactions among themselves. In the new 
sector we take the lightest stable particle to have no 
direct standard model coupling. The particle may still 
annihilate into standard model particles via lighter, un- 
stable exotic intermediary particles. If the intermediary 
is sufficiently light, an excess of galactic electrons and 
positrons is generated while baryonic production is kine- 
matically disallowed. 

The model has several other novel features. Since the 
dark and standard model sectors decouple early on, the 
temperature of the dark matter is generically different 
from the photon temperature at the time the relic abun- 
dance is formed. This modifies the allowed parameter 
space of relic masses and couplings. We also note that the 
dark-matter nucleon cross-section is much smaller than 
in many scenarios, which makes direct detection a chal- 
lenging prospect. 

In the following section we present a specific model 
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which we regard as a minimal existence proof of the 
feasabihty of reproducing the PAMELA positron excess 
from a dark sector containing a heavy stable particle and 
a light unstable intermediary, without any Somerfield en- 
hancement. The mechanism described can be embedded 
in a number of extensions to the Standard Model, in- 
cluding hidden valleys 23! and unparticles [25- Next, 
we examine the early-universe behavior of the model and 
derive the dark matter density today in terms of model 
parameters. We then calculate the present-day electron 
and positrons fluxes that result from galactic dark mat- 
ter annihilation, and show that they are consistent with 
PAMELA and ATIC. This model provides a reasonable 
fit to the shape of the Fermi spectrum, but produces 
more electrons than observed without a significant mod- 
ification to the boost factor. Finally, we will conclude 
with a few comments on some distinct features of this 
scenario. 



II. MODEL SETUP 

Consider the Standard Model extended by a number of 
new fields that we will label X, Y and Zi, . . . , Z„ where 
n > 0. The fields named Zi are unrelated to the Standard 
Model Z boson. For simplicity we will take all new fields 
to be scalars, though the model is can be easily modified 
to include fermionic dark matter or intermediaries. We'll 
take the Lagrangian to have the form 

C D M^X^ + m^Y'^ + m^Zf 

i 

+XX^Y^ + aYee + £,nt- (1) 

Here A and a are coupling constants, which we will treat 
as free parameters. The final term contains interactions, 
possibly strong, among the new particles so that none of 
the Zs are stable. We will require M ^ m > 2m^, and 
for simplicity we will take rrii ^ M so that the ZiS decay 
during the very early universe and their dynamics do not 
enter in the discussion below. The only interaction with 
the standard model fields is the coupling between the Y 
and electrons, which causes the Ys to decay at a rate 
Tye+e- • Note that the coupling a, which we take to be 
very small, is not gauge invariant under the electroweak 
transformations but could easily arise in a low energy 
effective theory below the scale of electroweak symme- 
try breaking. We neglect for now presentation of a UV 
complete model, which would address, for instance, other 
possible couplings and decays of the Y boson. 

Under these conditions, the X particle is stable and a 
WIMP dark matter candidate. The key quantity in both 
the early and late universe is the annihilation rate of Xs 
to Ys. The total cross-section for this process is given by 



where y/s is the center-of-mass energy. At temperatures 
well below M, the energy of the incoming Xs is 



« 2M + 3Tx 



(3) 



where Tx is the temperature of the X particles. To low- 
est order, the low-energy cross-section becomes 
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III. EARLY UNIVERSE 

The evolution of this model in the early universe is 
punctuated by several events. Note that the tempera- 
tures of the hidden and standard sectors are generally 
unequal after they decouple, which takes place at some 
early time idccoupic- We will assume that the Zi particles 
are heavy and have annihilated or decayed between times 
^decouple and txY , the time when the annihilation of X 
particles into Y particles goes out of thermal equilibrium. 
The main role of the Z particles is to allow for the pos- 
sibilility that the temperature in the hidden sector Tx is 
significantly higher than the temperature in the visible 
sector T-y at time txY, which will turn out to enhance the 
relic abundance of X particles today. The relic density 
of X particle dark matter is determined by the decou- 
pling of X and Y at time txY, which we assume to be 
later than ^decouple- We assume the Y particles are out of 
equilibrium with the standard model until time iydocayj 
given by the lifetime for Y decays, at which time the SM 
and the Y particles come back into thermal equilibrium. 
We will use the sudden decay approximation in which the 
reequilibration between the hidden and visible sectors oc- 
curs instantaneously. To avoid affecting the predictions 
of nucleosynthesis, we require that all these events occur 
before the nucleosynthesis epoch. 

After X decouples from Y, the X momentum simply 
redshifts with expansion, while the number of X parti- 
cles per co-moving volume is approximately conserved. 
Similarly, the entropy per co-moving volume is nearly 
conserved, except at time iYdocay It is thus convenient 
to compute the ratio of X particle density to entropy 
density. Once the increase in entropy density during the 
Y decay epoch is computed, we can determine with the 
dark matter to photon ratio today. 

We will begin with an approximate computation of the 
observed present-day abundance of dark matter in terms 
of the relic X energy density at txv- After txY, the 
dark matter energy density evolves as where R is 

the cosmological scale factor. Note that X particles are 
nonrelativistic at decoupling, and we take the mass of X 
to be constant. Assuming adiabatic evolution (no out of 
equilibrium processes) between txv and tYdocay, we can 
use entropy conservation to relate the photon tempera- 
ture T^,xY at txY to the photon temperature r-y^<Ydccay 
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just before tydc 



7, <Ydccay 
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where Rxy / ^i^Ydccay is the ratio of the scale factor at 
tydccay to the scalc factor at txv, g*,XY is the effec- 
tive number of relativistic species in equihbriuni with the 
photon at time txv and (?*.Ydccay is the effective number 
of relativistic species in equilibrium with the photon at 
time tydccay At time iydocay, using the sudden decay 
approximation, the photons get reheated to a tempera- 
ture T^,>Ydocay We define a parameter s, the fractional 
change in the entropy in the visible sector at ivdecay as 



/ -^7,>Ydccay \ 
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Note that s in general satisfies the inequality 
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where Thid.Ydccay denotes the temperature the hidden 
sector. The inequality is nearly saturated when the Y 
bosons are relativistic when they decay. 

We can now find the ratio of pxitxv), the dark matter 
density at txY to px{to), the dark matter density today. 
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This calculation differs from a standard calculation of the 
WIMP abundance today by the factors of s and 
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Note that we can use eqjTjto find an upper bound on the 
relic density compared with a standard relic abundance 
calculation in which all particles except the WIMP are in 
thermal equilbrium calculation. We define F such that 
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Saturating the inequality requires that the temperature 
in the hidden sector is much larger than the temperature 
of the visible sector, and that the Y particles decay while 



relativistic. If we take the photon temperature at the 
time of XY decouphng to be between 5 and 80 GeV, then 
g*,XY only includes standard model particles lighter than 
the weak bosons, and we may take g^,,xY = 86|. Assum- 
ing the Y particles decay when the photon temperature 
is between 10 and 100 MeV gives g*, Ydccay = Hf (includ- 
ing the contribution of the Y particles) . An approximate 
upper bound on F is therefore 



F < 46.6 . 

The X energy density at txY is 

Px,XY = Mxnx 

where 
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This may be related to the dark matter energy density 
px,Q today as 
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The current temperature of the photons, T^.o ~ 2.35 x 
10^ ^■^ GcV, and the total present dark matter abundance 
is fix ~ O.llft.^^, where h is the Hubble constant in units 
of 100(km/s)/Mpc. Taking the favored value h sa 0.70, 
the energy density of dark matter today is px,o ~ 5 x 
10~^GeVcm~'^. For a given value of he model dependent 
factor F we can relate M and Thid,j>fy- Note that the 
relation we differs from a standard decoupling calculation 
by the factor of F. 

We plot the ratio of the mass and temperature Thid,XY 
of X at freezout in figure [l] for a number of values of F, 
assuming the maximal value of g<t^xY = 86 j. 
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FIG. 1: The ratio of dark matter mass M to its temperature 
at decoupling Thid.xr for F = AQ.Q/i where i ={1,2,4} from 
top to bottom. 
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IV. DARK MATTER IN THE GALAXY 

After structure formation occurs, dark matter is gravi- 
tationally pulled into a structures around galaxies called 
halos. The density of the halo grows toward the center 
of the galaxy, and values are well above the mean value 
px- The annihilation rate, 



r — nav, 



(15) 



scales with the square of the number density n and so 
significantly increases in these regions. For a non rel- 
ativistic v, tru is approximately independent of w. In 
this section we will calculate the resulting electron and 
positron spectrums and compare to observations made 
by PAMELA and ATIC. Note that for a thermal relic, 
av is approximately determined by requiring that 
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In a standard WIMP annihilation calculation where all 
other particles are in thermal equilibrium with the pho- 
tons at decoupling the observed dark matter abundance 
today is obtained for 
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In our scenario this is modified by a factor of 
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Note that F is largest when T^.xy can be neglected rela- 

five to Thid.xy. The factor F j gl l^^-s a-n upper bound 
of about 5 so the upper bound on the annihilation cross 
section is approximately 



av < 1.5 X 10" 



-25 



cm 
sec 



(19) 



A larger annihilation cross section would imply an upper 
bound on the temperature of the hidden sector which is 
below Thid.xF together with some nonthermal produc- 
tion process for X particles, such as direct production in 
inflaton decay, or the mechanism in ref. |12| . 

For the sake of concreteness we will focus on the 
positron spectrum. The relevant process for positron pro- 
duction is XX — )■ YY — )■ 2e~^e~. We fix the branching 
fraction for decays of the intermediary Y particles to elec- 
tron positron pairs to unity. The flux at a given point in 
the galaxy are found by solving the diffusion equation of 
the positron density. 



(20) 



where K and b are energy-dependent diffusion and energy 
loss coefficients, and Q is the positron source term. The 



procedure to obtain the electron spectrum is analogous. 
A set of solutions to the diffusion equation in the case 
of a static distribution with disappearing flux outside a 
cylindrical region are well-known in the literature. A 
concise summary appears in |25| . For points near the 
solar system, the solution for the positron spectrum is 



Ve+ 1 (PeV 



M 



(21) 

In this expression _B is a dimensionless "boost" factor 
that accounts for sub-halo structure, v^^ is the positron 
velocity, b{E) is the energy loss coefficient, pq is the dark 
matter density near the solar system, and / is the halo 
function defined in terms of the diffusion length Xd ■ f{E) 
is the injection spectrum, given by. 
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where dNg+ / dE^+ for our model will be calculated below. 
Recall that to leading order the cross-section for 

XX — YY has temperature dependence a oc ^ . The 

temperature cancels with the temperature dependence of 

1 

the velocity, vx oc T^, leaving the product av tempera- 
ture independent to leading order. 

We now calculate the positron injection spectrum 
dN/dE created by X annihilation. In the galactic frame, 
the total energy of the colliding X particles is y/s = 2M, 
so the ejected relativistic Y particles have the Lorentz 
factor 
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In the center of mass frame of Y, it decays to 
e+e^ isotropically and monochromatically. Boosting the 
positrons back to the galactic frame yields an energy of 



E, 



M 
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(24) 



where 9 is the angle between the Y spatial momentum 
and positron spatial momentum. Since the emission of 
positrons in the frame of Y is isotropic, the above relation 
is readily converted to the injection spectrum by 



dN^+ _ dcosO dN^+ 
dE^+ dEf.+ dcosO' 

where dN^+fd cos 6 = 1/2. The result is 

dNf.+ m 



dEf,+ 2\/M^ — m'^Vm? — 4ml 



(25) 



(26) 



which is shown in figure [2j 

There are two sets of parameters needed for a calcu- 
lation of the observed positron and electron spectrums. 
The first are the parameters of the model we have de- 
scribed above. Assuming the dark matter consists of 



5 



0.0014 
0.0012 
0.0010 
0.0008 
0.0006 
0.0004 
0.0002 



200 



400 



600 



800 



1000 



FIG. 2: The energy spectrum dN/dE of the electrons as a 
function of the lab frame energy for M=800GeV, m=200MeV 



X particles, the relevant parameters are M, m, and A. 
The second set describe the mass distribution of the dark 
matter halo, and the diffusion of positrons and electrons. 
The parameters of the second set which are not known 
with precision are B and the factors that enter the diffu- 
sion length, Kq and 6. For the latter set, we will study 
the exemplar models from recent dark matter literature, 
which are outlined in [25] and [26]. We considered three 
standard halo profiles, NFW, Moore and IsoT, and three 
positron diffusion models, denoted "min," "med," and 
"max" . 

In practice, there is little dependence on the choice 
amongst the profile models, and here we will only present 
results from the NFW profile model. In addition the 
"med" and "max" propagation models yield very similar 
spectrums at Earth. The cross-section, and hence the 
flux, has very little dependence on the value of the Y 
mass m provided that it satisfied the kinematic bound 
on electron production. We will fix its value at 200 MeV. 

The values of B and A are degenerate in generating 
dark matter annihilation fluxes. The flux $ depends only 
on an overall scaling of BX^. We will require that the 
boost factor is within the bounds usually discussed in 
the literature, less than 20. Some studies suggest that 
even values of B above 2-3 are disfavored [17] ■ We will 
also require that A is sufficiently small that the F-electron 
coupling remains perturbative, which is true up to values 
of A ^ 10. We will find that we can generate spectra well 
above background levels without approaching either of 
these bounds. 

We fit our model to the recent PAMELA and ATIC 
observations, both individually and simultaneously. The 
PAMELA data which has been publicly discussed so 
far was presented as the ratio of the total number of 
positrons to the total number of positrons and electrons. 
This ratio departs background levels at 10 GeV, and rises 
to a value of 10% at 60 GeV. We fit the required value of 
BX^ required to generate such a spectrum as a function 
of M. The result is presented in figure [3] 

Here there is some slight tension between the param- 
eter B, X and the cross-section required to generate 



the relic abundance. Smaller B values require larger 
A, however this pushes us against the bound on av re- 
quired to generate the present abundance from a ther- 
mal relic given in eq. (19). For an 800 GeV X par- 



ticle and the "med" propagation model, the fit pro- 
duces B = 4.6 for A = 4. However, these values yield 
av — 1.4 X 10^^"* cm'^/s. We hit our bound of B = 20 for 
A — 1.9, which produces av = 3.3 x 10~^^ cm'^/s. This 
tension is significantly relaxed for smaller masses M. 




M (GeV) 

100 130 2UU JUU 3UU /OO lUOO 

FIG. 3; Required values of to generate the positron ex- 
cess described in the text. The top dashed curve is the "min" 
model, and the solid curve is the "med" models. 

Spectra for two values of M are shown for the lower 
half of the PAMELA energy range in figure |4] where the 
B factor has chosen to be optimal. The curves for two 
values of M are nearly degenerate below 60 GeV. Lower 
values of M show distinct spectral shapes that would be 
easily distinguished in the full PAMELA data release. 
Higher values of M have nearly identical shapes in the 
PAMELA range. 

Positron Ratio 
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0.05 



Energy (GeV) 



FIG. 4: Positron excess below 120GeV. The dash-dot curve 
is M = 100 GeV. The dashed curve is M = 800 GeV. The 
bottom dotted line is the background level. 

We now turn our attention to the ATIC results. The 
collaboration presented a spectrum in terms of total elec- 
tron flux that exhibits an excess well above their expected 
background level. The excess begins at 300 GeV and 
peaks at about 500 GeV. Their data at 900 GeV and 
above show no evidence for an excess. At the peak, the 
total number of electrons is around twice the background 
level. 
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We fit the ATIC data in terms of tlie ratio of the total 
number of electrons to the expected background num- 
ber. We find that the excess is well fit for values of 
M ^ 700 GeV. Remarkably, the best-fit boost values are 
quite similar to those required by the PAMELA data 
for masses in this range. The ATIC spectrum, using 
M = 800 GeV and the best-fit PAMELA boost factor, is 
shown in figure [5] 

N{£'')/Background 




N(e )/Background 




[ — , , 1 — , — , — ^ 1 , , 1 — , — , — J Energy{GeV) 

50 100 200 500 

FIG. 6: The ratio of the Fermi signal to background along 
with the spectrum for the dark matter model with parameters 
described in the text. Solid error bars indicate statistical 
errors, and dotted lines indicate the sum of statistical and 
systematic errors. 



' — ''''''' ' ■ ' ■ ' — ■ — ^ — ' — ■ — Energy (GeV) 

150 200 300 500 700 1000 

FIG. 5: The ratio of the ATIC signal to background along 
with the model fit for the parameters described in the text. 

We also fit the Fermi electron data fT' to our model. To 
estimate the background, we make a power law fit to the 
Fermi points below 100 GeV, below the bump feature. 
We then rescale the spectrum to this background, as we 
had for ATIC, and find the best-fit balues for M and 
boost. We found reasonable results for M = 700 GeV 
and boost = 0.4. The resulting spectrum is shown 
in figure [6j The spectrum generated by our model is 
softer than the Fermi spectrum when only statistical er- 
rors are considered, however the agreement greatly im- 
proves when systematic errors are included. Note that 
there is tension between the parameters which give a 
best fit to Fermi data and those that give the best fit 
to PAMELA / ATIC data in figure [s] Such tension is 
common to all models of direct annihilation of sub-TeV 
dark matter into messenger particles which decay into 
e+e" pairs, and can be reduced or eliminated by assum- 
ing a larger uncertainty in the Fermi energy resolution 



V. CONCLUSION 

Despite the compelling case for dark matter, we have 
few clues as to its nature, so it is important to consider 
a variety of possibilities with distinctive signatures. A 
nonminimal dark sector, consisting of more than one new 
particle interacting among themselves relatively strongly, 
occurs in a number of interesting scenarios, such as hid- 
den valley models [23] and unparticle models [23] . In such 
scenarios, a stable WIMP may annihilate into lighter 
unstable particles in the dark sector. The cross sec- 
tion for such annihilation can be larger than in typical 
WIMP models. The cross section for direct detection 
of dark matter would be very small, while decay of the 



lighter unstable particles into standard model particles 
offers new possibilities for indirect detection. Such en- 
hanced indirect detection signals are particularly inter- 
esting in light of the recent PAMELA and ATIC reports, 
as well as earlier reported excesses in the cosmic ray spec- 
trum. In this paper we have considered the cosmology 
of a simple example, where the light intermediary is a 
scalar particle, which may produce an excess of high en- 
ergy positrons and electrons in cosmic rays. We have 
shown that the model may easily explain the PAMELA 
and ATIC data. If the scalar coupling to fermions is 
proportional to fermion mass, as would be expected if 
this coupling arises from mixing with the Higgs, electron 
decay dominance would imply that the scalar is lighter 
than about 200 McV. Decays to 2 photons mediated by 
loops of heavier fermions would also be expected, giv- 
ing a gamma ray signal, possibly visible in EGRET and 
FERMI gamma ray observatory data [Ml ISO]- Such a 
gamma ray spectrum would distinguish our model from 
those with a light vector intermediary [T5til7j . as a vec- 
tor particle cannot decay to 2 photons, and decays to 
3 photons would be typically expected to have a small 
branching fraction. In addition, a scalar would not decay 
to neutrinos, so indirect detection would not be possible 
in neutrino observatories |31j . In the minimal version of 
the model we have presented, the dark sector would not 
show up in collider experiments, although it would be 
simple to extend the model so that additional new parti- 
cles with stronger standard model interactions could be 
produced at the LHC whose decay products include X 
or Y particles. 
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